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Abstract 

We summarize our results for hadronic contributions to the anomalous magnetic moment of 
the muon (a^), the one from hadronic vacuum-polarisation (HVP) and the light-by-light scattering 
contribution (LBL), obtained from the Dyson-Schwinger equations (DSE's) of QCD. In the case of 
pi I HVP we find good agreement with model independent determinations from dispersion relations for 

a^^^ as well as for the Adler function with deviations well below the ten percent level. From this 
D we conclude that the DSE approach should be capable of describing a]^^^ with similar accuracy. 



We also present results for LBL using a resonance expansion of the quark anti-quark T-matrix. 

^,1 



Our preliminary value is aj;^^ = (217 ± 91) x 10"^^ 



> 

o> 1 Introduction 

On 

• In the search for new physics beyond the standard model the anomalous magnetic moment of the muon 
^ (a^) is one of the most interesting observables. Compared to the corresponding electron anomaly (ae) 
^ it is more sensitive to contributions from high lying scales. These include the weak interactions, QCD 
and potential new physics [T]. Especially the contributions from soft QCD desire highest attention 
because, due to their non-perturbative nature and the resulting technical complications, they dominate 
^ the theoretical standard model (SM) prediction. 

5^ The efforts of the E821 experiment at Brookhaven National Lab [21 [3] as well as theoretical efforts of 
more than a decade [1] culminated in a determination of down to a level where significant deviations 
have been found 

Experiment: 116 592 089(63) x 10"" , (1) 
Theory: 116 591 828(49) x 10"" , (2) 

where the theoretical number is taken from Ref. 0. Comparing theory and experiment the deviation 
amounts to a^^P — a^^*^" = 261(80) which corresponds to a 3.2a effect. In order to confirm this result 
the uncertainties have to be reduced further. 

There are two hadronic contributions that dominate the SM uncertainty. There is the hadronic 
vacuum polarisation contribution (HVP) which gives rise to the leading hadronic contribution as well 
as the leading SM uncertainty contribution [S] 

^HVP,DR ^ g 949.1(42.7) X 10"". (3) 
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(a) (b) 

Figure 1: The two classifications of corrections to tlie plioton-muon vertex function: (a) hadronic 
vacuum polarization contribution to a^. The vertex is dressed by the vacuum polarization tensor n^j,; 
(b) the hadronic light-by-light scattering contribution to a^. 



The relevant diagram, involving the hadronic one-particle irreducible (IPI) photon self-energy II^j, is 



shown in figure 1(a) The next-to leading uncertainty contribution comes from the light-by-light (LBL) 
scattering contribution that is shown in fig. l(b)[ Estimates from the viewpoint of effective field theory 
(EFT) from different approaches were recently combined into a single number [6] 



,LBL 



105(26) X 10 



-11 



(4) 



The uncertainty given here is rather small compared to most estimates. In fact our results indicate that 
this error may be far too optimistic. 

Our strategy to determine these quantities is the following. We work with the Dyson- Schwinger and 
Bethe-Salpeter equations (DSE/BSE) of QCD [3,18]. With these we calculate the HVP contribution to 

where we use a parameter set, among others, that is completely fixed by meson phenomenology. The 
HVP contribution can be compared to essentially model independent result from dispersion relations 
[H] such that the calculation serves as a non-trivial cross check of our methods. Afterwards we approach 
the LBL contribution using exactly the same truncation such that we have reasons to believe that we 
can ultimately reach a similar precision as in the case of HVP. 

This proceedings contribution is organized as follows. First of all we summarise the employed 
truncation in section |2l The HVP contribution will be discussed in section |3] and LBL in sec. HI 
Afterwards we discuss our results for both of these contributions in section |5] and conclude. 



2 Calculational scheme 

We work in rainbow-ladder truncation of QCD using the Maris- Tandy model of the quark-gluon inter- 
action [TU]. The central object in this approach is the quark DSE 

S-\p) = Z, So' + l,S{ph, T,,{k) G(P), (5) 

<? 

where S is the full quark propagator, 5*0 the corresponding bare quantity and Z2 is the quark wave- 
function renormalisation. T^y{k) is the transverse projector and G{k'^) is the effective gluon dressing. 
This function is modelled in the present approach in a way such that chiral symmetry breaking occurs 
while the axial- vector Ward-Takahashi identity (AXWTI), the U{1) vector- WTI of QED and resummed 
one-loop perturbation theory are respected lilOj. Consistent with the quark DSE in (|5| is the meson 
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BSE 

[TUP, k) = ~Zl^ J [S(g+)r(P, q)Siq.)UKtu,rsik - q), (6) 
g 

where P is the meson momentum, k the relative quark momentum and q± = q ± P/2. The interaction 
kernel is defined as 

Krs,tu{k) = Gik')T,,{k)[^,]J^,]^^. (7) 

The latter two equations are intimately related by chiral symmetry, to give a dynamical breaking in 
accordance with Goldstone's theorem [H]. In addition we need the quark-photon vertex defined via the 
inhomogeneous BSE 

[T,UP, k) = Z,^, [S{q+)^,{P. q)S{q-)]utKtuAk - q), (8) 

which is the key to any calculation of electromagnetic properties of hadrons. The vertex features 
a vector-meson bound-state for time-like momenta P such that vector-meson dominance (VMD) is 
dynamically included. This gives e.g. important contributions to the pion charge radius which can be 
nicely described in the present approach |T2]. 



3 Hadronic vacuum polarisation (HVP) 

Here we present our results briefiy, more details can be found in Ref. pL3j. The central object for the 
HVP contribution is the hadronic tensor 

U^p) = -Z^e'J Tr[5(g+)r^(p,g)5(g_)7.] , (9) 

which corresponds to the IPI hadronic photon self-energy and involves the non-perturbative quark 
propagator ([s]) and the self-consistent quark-photon vertex (|8|. The tensor H^,^ is transverse due to its 
WTI 



n,.(p)=(5,.-^)/H(/) , (10) 

which serves as a definition of the scalar function H(p^). We use the renormalisation condition Hr(p^) = 
H(p^) — H(0) which corresponds to the usual physical QED on-shell scheme giving rise to e(0) = Cphysicai- 
Another quantity that is interesting in the present context is the Adler function 

z,(„ = . (n) 

The results for the Adler function from DSE's including five quark fiavours is shown in fig. [2] and 
compared to a result from dispersion relations [9J. There is quite reasonable agreement at all momentum 
scales. Especially the deeply non-perturbative behaviour below 1 GeV is nicely reproduced. It is this 
regime, set by the muon mass, where the contribution to saturates. 
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Figure 2: The Adler function obtained from DSE's using our parameter set II compared to a model 
independent result from dispersion relations. 



[MeV] 




rris 




rriK 






set I 


3.7 


85 


138 


495 


740 


1080 


set II 


11 


72 


240 


477 


770 


1020 



Table 1: Two choices for the light bare quark masses at = (19 GeV)^ and the resulting meson masses 
(in MeV) in the pseudoscalar and vector meson sector. For the heavy quarks we always take rric = 827 
MeV and = 3680 MeV which lead to good results for charmonia and bottomonia in the pseudoscalar 
and vector channel. 



The contribution to the muon g — 2 can be obtained via 



HYP 



a 



dx (1 — x) 



-e'nR 



1-x ^ 



(12) 



where a is the fine-structure constant. In particular we use two different parameter sets. The standard 
parameter set where the u, d and s quark masses are fixed to the pseudo-scalar meson sector and 
another one where instead the vector-meson sector is used. This is summarized in table [U The c and 
b quark mass functions are fixed to the charmonium and bottomonium vector meson states in all cases 
|15j . With these two parameter sets we obtain the following results 



,HVPJ 



7440 X 10 



-11 



HVPJI 



6760 X 10 



-11 



(13) 



Comparing to the model independent result ^ we see that the standard set I deviates about 7 percent. 
We see the reason for this in the p mass which is about four percent too light, see tab. [T} The result 
with the physical p mass (set II) is indeed closer to the dispersion relation result. Taking the idea of 
changing the p mass a step further, we calculate the two fiavour contribution function 
of the vector meson mass. Both are functions of the quark masses m^/^. In addition we calculate 
where the strange contribution is just an additive constant since rris remains fixed. The 



results are shown in fig. [3] where we compare to Nf = 2 results from the ETMC collaboration (TB] (blue 
data) and Nf = 2 + 1 data (red) from the RBC-UKQCD collaboration [T7|. The ordinate shows in 
units of [10^"]. Our curves agree with both data within error bars which we take as a hint that the 
DSE/BSE approach in the present truncation captures the relevant degrees of freedom. For HVP this 
seems to be more than anything else the vector meson as would be expected from VMD estimates [18] . 
For a detailed discussion see Ref. [131 . 
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Figure 3: The Nf = 2 and 



1 flavour contribution to units of [10 ^^1 function of 



tlie vector meson mass. For tlie latter case the s contribution is kept constant. The two DSE curves 
are compared to lattice results for the two- and three flavour case respectively. The blue data is iV/ = 2 
data from ETMC [IS] and the red data is 2 + 1 from RBC-UKQCD [H]. 

4 Hadronic light- by- light scattering (LBL) 

In the present section we discuss the LBL contribution. Within the framework of DSE's the hadronic 
four-point function, that is the essential ingredient here, has a description that is consistent with the one 
for HVP shown earlier We presented this representation in Refs. [191 EQ] where also more details 
can be found. There we also elaborate on the resonance expansion of the quark anti-quark T-matrix 
that is used for the results presented here. To this end we arrive at an approximate description of the 
full four-point function 





(14) 



that consists of the non-perturbatively dressed quark loop diagram (QL) as well as a pseudo-scalar 
meson-exchange contribution that takes into account the 7r°, rj and t]' mesons. This picture is very 
similar to the one obtained in hadronic models and EFT p31|2T| [221 [2H| [2it [25| |26| . 
In order to obtain the contribution to a^^^ from the four-point function we define 



ieT 




(15) 



from which the anomaly can be obtained via [27] 

1 



48m, 



-Tr 



{if" + m^)[7^,7p](i/' + m^)V 



ap 



(16) 



k=0 



where P is the muon momentum. Here 5* and D^j^ are perturbative muon and photon propagators and 
the definition Il(p)fj,uai3 = dp^^uap has been used with the hadronic four-point function li^^ap- 
For the pseudo-scalar (PS) meson pole contribution we need the PS77 form factor 



kll^\h,k2) = 2e^Nj Tr[tQeT.ik2,Pl2)SFiP2)T''''ip23,P)SFiP3)tQeT^ih,Psi)SFip2 



(17) 
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that is defined as a non-perturbative quark triangle that involves the quark ([s]), the quark-photon 
vertex (js) and the meson amplitude (6) called f^^ here. The meson momentum is P, ki/2 are the 
photon momenta, pi the momenta of the quarks and pij = {pi + Pj) /2. From the form factor together 
with a bare meson propagator we obtain the resonant part that is shown in Eq. (14). Details concerning 



the flavour content of the meson as well as the necessary meson off-shell prescription can be found in 
|20j . Once the form factor is known the contribution to g — 2 can be obtained along the lines explained 
in [21]. Our result for the PS meson exchange contribution (vr^, t], t]') is 



^LBL;PS 



(80.7 ± 12.0) X 10 



-11 



where the error is dominantly an estimate of the systematic model uncertainty. 

For the QL contribution we take the full quark propagator ([s]) together with the Ball-Chiu (BC) 
vertex construction [28] 



where the symbols 



F{ki) + F{e_] 



F{kl) - F{kl 
kl-k^ 



(19) 



(20) 



have been used. This substructure of the fully self-consistent vertex (|8j) is dictated by the vector WTI. 
Defining the four-point function from the quark-loop as in Eq. (14), taking the derivative and using 



Eqs. (15) and (16) we obtain 



a 



a 



a 



LBL;quarkloop (bare vertex) 
LBL;quarkloop (IBC) 
LBL;quarkloop (BC) 



( 61 ±2) X 10-11 
(107 ± 2) X 10-11 
(176 ±4) X 10-11, 



(21) 



where the first result uses a bare quark-photon vertex and IBC only has the 7^ part of the dressed 



vertex, see (19). In this calculation we included the flavours m, c?, s and c. It can clearly be seen that the 



vertex dressing causes quite some enhancement especially when all three structures of the BC vertex 
are used (third case). The error given here is numerical. Clearly, this result is preliminary, since the 
transverse structure of the vertex including important contributions from vector mesons is missing. We 
have estimated these contributions in [20] , however a full calculation is absolutely mandatory and well 
under way. 



5 Discussion 

We presented results for the HVP as well as for the LBL contribution to the muon g — 2 obtained 



within the framework of DSE's. We saw that in the case of HVP our results for a^^^ (13) reproduce 
model independent dispersion relations on the less-than-ten-percent level. We see no principal reason 
why this should not be the case also for a full LBL calculation. Indeed, our result for the pseudo- 
scalar meson-exchange contribution to LBL ( [Is] ) is in the ballpark of the results obtained within other 
approaches [TH EU |22l |23l [211 EHl EE] (see [4j for an overview). For the quark loop contribution 
to LBL we take our BC result from ( [2l] ) and guesstimate the missing contributions from the above 
mentioned transverse components of the quark-photon vertex, see [2^ for details. Thus we arrive at 
the value gj|^^L,quarkioop (Bc+transverse) _ ^^gg _|_ ^g-^ ^ iQ-n, whcrc the large error takes into account the 

uncertainties due to the estimate. Puting all contributions together we obtain 

^LBL ^ (217 ± 91) X 10-11 , (22) 
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As mentioned already our result for LBL hints towards a larger contribution and thus to a smaller 
deviation between theory and experiment as compared to Eqs. [2]). Besides implementing the full 
quark-photon vertex inside the quark-loop, another important next step is to overcome the resonance 
approximation of LBL. In general, we believe to have shown that a full calculation from the DSE 
approach can be expected to mark a clear step forward for the case of the hadronic light-by-light 
scattering contribution to the muon g — 2. 
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